\(\int \frac {\sqrt {a+b x}}{\sqrt {c+d x} \sqrt {e+f x} \sqrt {g+h x}} \, dx\) [108]

   Optimal result
   Rubi [A] (verified)
   Mathematica [B] (verified)
   Maple [B] (verified)
   Fricas [F(-1)]
   Sympy [F]
   Maxima [F]
   Giac [F]
   Mupad [F(-1)]

Optimal result

Integrand size = 37, antiderivative size = 228 \[ \int \frac {\sqrt {a+b x}}{\sqrt {c+d x} \sqrt {e+f x} \sqrt {g+h x}} \, dx=\frac {2 \sqrt {-d g+c h} (a+b x) \sqrt {\frac {(b g-a h) (c+d x)}{(d g-c h) (a+b x)}} \sqrt {\frac {(b g-a h) (e+f x)}{(f g-e h) (a+b x)}} \operatorname {EllipticPi}\left (-\frac {b (d g-c h)}{(b c-a d) h},\arcsin \left (\frac {\sqrt {b c-a d} \sqrt {g+h x}}{\sqrt {-d g+c h} \sqrt {a+b x}}\right ),\frac {(b e-a f) (d g-c h)}{(b c-a d) (f g-e h)}\right )}{\sqrt {b c-a d} h \sqrt {c+d x} \sqrt {e+f x}} \]

[Out]

2*(b*x+a)*EllipticPi((-a*d+b*c)^(1/2)*(h*x+g)^(1/2)/(c*h-d*g)^(1/2)/(b*x+a)^(1/2),-b*(-c*h+d*g)/(-a*d+b*c)/h,(
(-a*f+b*e)*(-c*h+d*g)/(-a*d+b*c)/(-e*h+f*g))^(1/2))*(c*h-d*g)^(1/2)*((-a*h+b*g)*(d*x+c)/(-c*h+d*g)/(b*x+a))^(1
/2)*((-a*h+b*g)*(f*x+e)/(-e*h+f*g)/(b*x+a))^(1/2)/h/(-a*d+b*c)^(1/2)/(d*x+c)^(1/2)/(f*x+e)^(1/2)

Rubi [A] (verified)

Time = 0.10 (sec) , antiderivative size = 228, normalized size of antiderivative = 1.00, number of steps used = 2, number of rules used = 2, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.054, Rules used = {171, 551} \[ \int \frac {\sqrt {a+b x}}{\sqrt {c+d x} \sqrt {e+f x} \sqrt {g+h x}} \, dx=\frac {2 (a+b x) \sqrt {c h-d g} \sqrt {\frac {(c+d x) (b g-a h)}{(a+b x) (d g-c h)}} \sqrt {\frac {(e+f x) (b g-a h)}{(a+b x) (f g-e h)}} \operatorname {EllipticPi}\left (-\frac {b (d g-c h)}{(b c-a d) h},\arcsin \left (\frac {\sqrt {b c-a d} \sqrt {g+h x}}{\sqrt {c h-d g} \sqrt {a+b x}}\right ),\frac {(b e-a f) (d g-c h)}{(b c-a d) (f g-e h)}\right )}{h \sqrt {c+d x} \sqrt {e+f x} \sqrt {b c-a d}} \]

[In]

Int[Sqrt[a + b*x]/(Sqrt[c + d*x]*Sqrt[e + f*x]*Sqrt[g + h*x]),x]

[Out]

(2*Sqrt[-(d*g) + c*h]*(a + b*x)*Sqrt[((b*g - a*h)*(c + d*x))/((d*g - c*h)*(a + b*x))]*Sqrt[((b*g - a*h)*(e + f
*x))/((f*g - e*h)*(a + b*x))]*EllipticPi[-((b*(d*g - c*h))/((b*c - a*d)*h)), ArcSin[(Sqrt[b*c - a*d]*Sqrt[g +
h*x])/(Sqrt[-(d*g) + c*h]*Sqrt[a + b*x])], ((b*e - a*f)*(d*g - c*h))/((b*c - a*d)*(f*g - e*h))])/(Sqrt[b*c - a
*d]*h*Sqrt[c + d*x]*Sqrt[e + f*x])

Rule 171

Int[Sqrt[(a_.) + (b_.)*(x_)]/(Sqrt[(c_.) + (d_.)*(x_)]*Sqrt[(e_.) + (f_.)*(x_)]*Sqrt[(g_.) + (h_.)*(x_)]), x_S
ymbol] :> Dist[2*(a + b*x)*Sqrt[(b*g - a*h)*((c + d*x)/((d*g - c*h)*(a + b*x)))]*(Sqrt[(b*g - a*h)*((e + f*x)/
((f*g - e*h)*(a + b*x)))]/(Sqrt[c + d*x]*Sqrt[e + f*x])), Subst[Int[1/((h - b*x^2)*Sqrt[1 + (b*c - a*d)*(x^2/(
d*g - c*h))]*Sqrt[1 + (b*e - a*f)*(x^2/(f*g - e*h))]), x], x, Sqrt[g + h*x]/Sqrt[a + b*x]], x] /; FreeQ[{a, b,
 c, d, e, f, g, h}, x]

Rule 551

Int[1/(((a_) + (b_.)*(x_)^2)*Sqrt[(c_) + (d_.)*(x_)^2]*Sqrt[(e_) + (f_.)*(x_)^2]), x_Symbol] :> Simp[(1/(a*Sqr
t[c]*Sqrt[e]*Rt[-d/c, 2]))*EllipticPi[b*(c/(a*d)), ArcSin[Rt[-d/c, 2]*x], c*(f/(d*e))], x] /; FreeQ[{a, b, c,
d, e, f}, x] &&  !GtQ[d/c, 0] && GtQ[c, 0] && GtQ[e, 0] &&  !( !GtQ[f/e, 0] && SimplerSqrtQ[-f/e, -d/c])

Rubi steps \begin{align*} \text {integral}& = \frac {\left (2 (a+b x) \sqrt {\frac {(b g-a h) (c+d x)}{(d g-c h) (a+b x)}} \sqrt {\frac {(b g-a h) (e+f x)}{(f g-e h) (a+b x)}}\right ) \text {Subst}\left (\int \frac {1}{\left (h-b x^2\right ) \sqrt {1+\frac {(b c-a d) x^2}{d g-c h}} \sqrt {1+\frac {(b e-a f) x^2}{f g-e h}}} \, dx,x,\frac {\sqrt {g+h x}}{\sqrt {a+b x}}\right )}{\sqrt {c+d x} \sqrt {e+f x}} \\ & = \frac {2 \sqrt {-d g+c h} (a+b x) \sqrt {\frac {(b g-a h) (c+d x)}{(d g-c h) (a+b x)}} \sqrt {\frac {(b g-a h) (e+f x)}{(f g-e h) (a+b x)}} \Pi \left (-\frac {b (d g-c h)}{(b c-a d) h};\sin ^{-1}\left (\frac {\sqrt {b c-a d} \sqrt {g+h x}}{\sqrt {-d g+c h} \sqrt {a+b x}}\right )|\frac {(b e-a f) (d g-c h)}{(b c-a d) (f g-e h)}\right )}{\sqrt {b c-a d} h \sqrt {c+d x} \sqrt {e+f x}} \\ \end{align*}

Mathematica [B] (verified)

Leaf count is larger than twice the leaf count of optimal. \(583\) vs. \(2(228)=456\).

Time = 31.71 (sec) , antiderivative size = 583, normalized size of antiderivative = 2.56 \[ \int \frac {\sqrt {a+b x}}{\sqrt {c+d x} \sqrt {e+f x} \sqrt {g+h x}} \, dx=-\frac {2 \sqrt {\frac {(d g-c h) (a+b x)}{(b g-a h) (c+d x)}} (c+d x)^{3/2} \left (\frac {a d \sqrt {\frac {(d g-c h) (e+f x)}{(f g-e h) (c+d x)}} (g+h x) \operatorname {EllipticF}\left (\arcsin \left (\sqrt {\frac {(-d e+c f) (g+h x)}{(f g-e h) (c+d x)}}\right ),\frac {(b c-a d) (-f g+e h)}{(d e-c f) (b g-a h)}\right )}{(d g-c h) (c+d x) \sqrt {\frac {(-d e+c f) (g+h x)}{(f g-e h) (c+d x)}}}+\frac {b c \sqrt {\frac {(d g-c h) (e+f x)}{(f g-e h) (c+d x)}} (g+h x) \operatorname {EllipticF}\left (\arcsin \left (\sqrt {\frac {(-d e+c f) (g+h x)}{(f g-e h) (c+d x)}}\right ),\frac {(b c-a d) (-f g+e h)}{(d e-c f) (b g-a h)}\right )}{(-d g+c h) (c+d x) \sqrt {\frac {(-d e+c f) (g+h x)}{(f g-e h) (c+d x)}}}+\frac {b (f g-e h) \sqrt {\frac {(-d e+c f) (d g-c h) (e+f x) (g+h x)}{(f g-e h)^2 (c+d x)^2}} \operatorname {EllipticPi}\left (\frac {d (-f g+e h)}{(d e-c f) h},\arcsin \left (\sqrt {\frac {(-d e+c f) (g+h x)}{(f g-e h) (c+d x)}}\right ),\frac {(b c-a d) (-f g+e h)}{(d e-c f) (b g-a h)}\right )}{(d e-c f) h}\right )}{d \sqrt {a+b x} \sqrt {e+f x} \sqrt {g+h x}} \]

[In]

Integrate[Sqrt[a + b*x]/(Sqrt[c + d*x]*Sqrt[e + f*x]*Sqrt[g + h*x]),x]

[Out]

(-2*Sqrt[((d*g - c*h)*(a + b*x))/((b*g - a*h)*(c + d*x))]*(c + d*x)^(3/2)*((a*d*Sqrt[((d*g - c*h)*(e + f*x))/(
(f*g - e*h)*(c + d*x))]*(g + h*x)*EllipticF[ArcSin[Sqrt[((-(d*e) + c*f)*(g + h*x))/((f*g - e*h)*(c + d*x))]],
((b*c - a*d)*(-(f*g) + e*h))/((d*e - c*f)*(b*g - a*h))])/((d*g - c*h)*(c + d*x)*Sqrt[((-(d*e) + c*f)*(g + h*x)
)/((f*g - e*h)*(c + d*x))]) + (b*c*Sqrt[((d*g - c*h)*(e + f*x))/((f*g - e*h)*(c + d*x))]*(g + h*x)*EllipticF[A
rcSin[Sqrt[((-(d*e) + c*f)*(g + h*x))/((f*g - e*h)*(c + d*x))]], ((b*c - a*d)*(-(f*g) + e*h))/((d*e - c*f)*(b*
g - a*h))])/((-(d*g) + c*h)*(c + d*x)*Sqrt[((-(d*e) + c*f)*(g + h*x))/((f*g - e*h)*(c + d*x))]) + (b*(f*g - e*
h)*Sqrt[((-(d*e) + c*f)*(d*g - c*h)*(e + f*x)*(g + h*x))/((f*g - e*h)^2*(c + d*x)^2)]*EllipticPi[(d*(-(f*g) +
e*h))/((d*e - c*f)*h), ArcSin[Sqrt[((-(d*e) + c*f)*(g + h*x))/((f*g - e*h)*(c + d*x))]], ((b*c - a*d)*(-(f*g)
+ e*h))/((d*e - c*f)*(b*g - a*h))])/((d*e - c*f)*h)))/(d*Sqrt[a + b*x]*Sqrt[e + f*x]*Sqrt[g + h*x])

Maple [B] (verified)

Leaf count of result is larger than twice the leaf count of optimal. \(847\) vs. \(2(209)=418\).

Time = 1.29 (sec) , antiderivative size = 848, normalized size of antiderivative = 3.72

method result size
elliptic \(\frac {\sqrt {\left (b x +a \right ) \left (d x +c \right ) \left (f x +e \right ) \left (h x +g \right )}\, \left (\frac {2 a \left (\frac {g}{h}-\frac {a}{b}\right ) \sqrt {\frac {\left (-\frac {g}{h}+\frac {c}{d}\right ) \left (x +\frac {a}{b}\right )}{\left (-\frac {g}{h}+\frac {a}{b}\right ) \left (x +\frac {c}{d}\right )}}\, \left (x +\frac {c}{d}\right )^{2} \sqrt {\frac {\left (-\frac {c}{d}+\frac {a}{b}\right ) \left (x +\frac {e}{f}\right )}{\left (-\frac {e}{f}+\frac {a}{b}\right ) \left (x +\frac {c}{d}\right )}}\, \sqrt {\frac {\left (-\frac {c}{d}+\frac {a}{b}\right ) \left (x +\frac {g}{h}\right )}{\left (-\frac {g}{h}+\frac {a}{b}\right ) \left (x +\frac {c}{d}\right )}}\, F\left (\sqrt {\frac {\left (-\frac {g}{h}+\frac {c}{d}\right ) \left (x +\frac {a}{b}\right )}{\left (-\frac {g}{h}+\frac {a}{b}\right ) \left (x +\frac {c}{d}\right )}}, \sqrt {\frac {\left (\frac {e}{f}-\frac {c}{d}\right ) \left (\frac {g}{h}-\frac {a}{b}\right )}{\left (-\frac {a}{b}+\frac {e}{f}\right ) \left (-\frac {c}{d}+\frac {g}{h}\right )}}\right )}{\left (-\frac {g}{h}+\frac {c}{d}\right ) \left (-\frac {c}{d}+\frac {a}{b}\right ) \sqrt {b d f h \left (x +\frac {a}{b}\right ) \left (x +\frac {c}{d}\right ) \left (x +\frac {e}{f}\right ) \left (x +\frac {g}{h}\right )}}+\frac {2 b \left (\frac {g}{h}-\frac {a}{b}\right ) \sqrt {\frac {\left (-\frac {g}{h}+\frac {c}{d}\right ) \left (x +\frac {a}{b}\right )}{\left (-\frac {g}{h}+\frac {a}{b}\right ) \left (x +\frac {c}{d}\right )}}\, \left (x +\frac {c}{d}\right )^{2} \sqrt {\frac {\left (-\frac {c}{d}+\frac {a}{b}\right ) \left (x +\frac {e}{f}\right )}{\left (-\frac {e}{f}+\frac {a}{b}\right ) \left (x +\frac {c}{d}\right )}}\, \sqrt {\frac {\left (-\frac {c}{d}+\frac {a}{b}\right ) \left (x +\frac {g}{h}\right )}{\left (-\frac {g}{h}+\frac {a}{b}\right ) \left (x +\frac {c}{d}\right )}}\, \left (-\frac {c F\left (\sqrt {\frac {\left (-\frac {g}{h}+\frac {c}{d}\right ) \left (x +\frac {a}{b}\right )}{\left (-\frac {g}{h}+\frac {a}{b}\right ) \left (x +\frac {c}{d}\right )}}, \sqrt {\frac {\left (\frac {e}{f}-\frac {c}{d}\right ) \left (\frac {g}{h}-\frac {a}{b}\right )}{\left (-\frac {a}{b}+\frac {e}{f}\right ) \left (-\frac {c}{d}+\frac {g}{h}\right )}}\right )}{d}+\left (\frac {c}{d}-\frac {a}{b}\right ) \Pi \left (\sqrt {\frac {\left (-\frac {g}{h}+\frac {c}{d}\right ) \left (x +\frac {a}{b}\right )}{\left (-\frac {g}{h}+\frac {a}{b}\right ) \left (x +\frac {c}{d}\right )}}, \frac {-\frac {g}{h}+\frac {a}{b}}{-\frac {g}{h}+\frac {c}{d}}, \sqrt {\frac {\left (\frac {e}{f}-\frac {c}{d}\right ) \left (\frac {g}{h}-\frac {a}{b}\right )}{\left (-\frac {a}{b}+\frac {e}{f}\right ) \left (-\frac {c}{d}+\frac {g}{h}\right )}}\right )\right )}{\left (-\frac {g}{h}+\frac {c}{d}\right ) \left (-\frac {c}{d}+\frac {a}{b}\right ) \sqrt {b d f h \left (x +\frac {a}{b}\right ) \left (x +\frac {c}{d}\right ) \left (x +\frac {e}{f}\right ) \left (x +\frac {g}{h}\right )}}\right )}{\sqrt {b x +a}\, \sqrt {d x +c}\, \sqrt {f x +e}\, \sqrt {h x +g}}\) \(848\)
default \(\text {Expression too large to display}\) \(1250\)

[In]

int((b*x+a)^(1/2)/(d*x+c)^(1/2)/(f*x+e)^(1/2)/(h*x+g)^(1/2),x,method=_RETURNVERBOSE)

[Out]

((b*x+a)*(d*x+c)*(f*x+e)*(h*x+g))^(1/2)/(b*x+a)^(1/2)/(d*x+c)^(1/2)/(f*x+e)^(1/2)/(h*x+g)^(1/2)*(2*a*(g/h-a/b)
*((-g/h+c/d)*(x+a/b)/(-g/h+a/b)/(x+c/d))^(1/2)*(x+c/d)^2*((-c/d+a/b)*(x+e/f)/(-e/f+a/b)/(x+c/d))^(1/2)*((-c/d+
a/b)*(x+g/h)/(-g/h+a/b)/(x+c/d))^(1/2)/(-g/h+c/d)/(-c/d+a/b)/(b*d*f*h*(x+a/b)*(x+c/d)*(x+e/f)*(x+g/h))^(1/2)*E
llipticF(((-g/h+c/d)*(x+a/b)/(-g/h+a/b)/(x+c/d))^(1/2),((e/f-c/d)*(g/h-a/b)/(-a/b+e/f)/(-c/d+g/h))^(1/2))+2*b*
(g/h-a/b)*((-g/h+c/d)*(x+a/b)/(-g/h+a/b)/(x+c/d))^(1/2)*(x+c/d)^2*((-c/d+a/b)*(x+e/f)/(-e/f+a/b)/(x+c/d))^(1/2
)*((-c/d+a/b)*(x+g/h)/(-g/h+a/b)/(x+c/d))^(1/2)/(-g/h+c/d)/(-c/d+a/b)/(b*d*f*h*(x+a/b)*(x+c/d)*(x+e/f)*(x+g/h)
)^(1/2)*(-c/d*EllipticF(((-g/h+c/d)*(x+a/b)/(-g/h+a/b)/(x+c/d))^(1/2),((e/f-c/d)*(g/h-a/b)/(-a/b+e/f)/(-c/d+g/
h))^(1/2))+(c/d-a/b)*EllipticPi(((-g/h+c/d)*(x+a/b)/(-g/h+a/b)/(x+c/d))^(1/2),(-g/h+a/b)/(-g/h+c/d),((e/f-c/d)
*(g/h-a/b)/(-a/b+e/f)/(-c/d+g/h))^(1/2))))

Fricas [F(-1)]

Timed out. \[ \int \frac {\sqrt {a+b x}}{\sqrt {c+d x} \sqrt {e+f x} \sqrt {g+h x}} \, dx=\text {Timed out} \]

[In]

integrate((b*x+a)^(1/2)/(d*x+c)^(1/2)/(f*x+e)^(1/2)/(h*x+g)^(1/2),x, algorithm="fricas")

[Out]

Timed out

Sympy [F]

\[ \int \frac {\sqrt {a+b x}}{\sqrt {c+d x} \sqrt {e+f x} \sqrt {g+h x}} \, dx=\int \frac {\sqrt {a + b x}}{\sqrt {c + d x} \sqrt {e + f x} \sqrt {g + h x}}\, dx \]

[In]

integrate((b*x+a)**(1/2)/(d*x+c)**(1/2)/(f*x+e)**(1/2)/(h*x+g)**(1/2),x)

[Out]

Integral(sqrt(a + b*x)/(sqrt(c + d*x)*sqrt(e + f*x)*sqrt(g + h*x)), x)

Maxima [F]

\[ \int \frac {\sqrt {a+b x}}{\sqrt {c+d x} \sqrt {e+f x} \sqrt {g+h x}} \, dx=\int { \frac {\sqrt {b x + a}}{\sqrt {d x + c} \sqrt {f x + e} \sqrt {h x + g}} \,d x } \]

[In]

integrate((b*x+a)^(1/2)/(d*x+c)^(1/2)/(f*x+e)^(1/2)/(h*x+g)^(1/2),x, algorithm="maxima")

[Out]

integrate(sqrt(b*x + a)/(sqrt(d*x + c)*sqrt(f*x + e)*sqrt(h*x + g)), x)

Giac [F]

\[ \int \frac {\sqrt {a+b x}}{\sqrt {c+d x} \sqrt {e+f x} \sqrt {g+h x}} \, dx=\int { \frac {\sqrt {b x + a}}{\sqrt {d x + c} \sqrt {f x + e} \sqrt {h x + g}} \,d x } \]

[In]

integrate((b*x+a)^(1/2)/(d*x+c)^(1/2)/(f*x+e)^(1/2)/(h*x+g)^(1/2),x, algorithm="giac")

[Out]

integrate(sqrt(b*x + a)/(sqrt(d*x + c)*sqrt(f*x + e)*sqrt(h*x + g)), x)

Mupad [F(-1)]

Timed out. \[ \int \frac {\sqrt {a+b x}}{\sqrt {c+d x} \sqrt {e+f x} \sqrt {g+h x}} \, dx=\int \frac {\sqrt {a+b\,x}}{\sqrt {e+f\,x}\,\sqrt {g+h\,x}\,\sqrt {c+d\,x}} \,d x \]

[In]

int((a + b*x)^(1/2)/((e + f*x)^(1/2)*(g + h*x)^(1/2)*(c + d*x)^(1/2)),x)

[Out]

int((a + b*x)^(1/2)/((e + f*x)^(1/2)*(g + h*x)^(1/2)*(c + d*x)^(1/2)), x)